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1. Thermodynamics of correlations in quark matter 

While numerous works have recently studied the thermodynamics of quark matter on the 
mean-field level including the effects of the medium dependence of the order parameters, not so 
much is known beyond the mean field, about hadronic correlations and their backreaction to the 
structure of the model QCD phase diagram and its thermodynamics. Here we will elaborate on the 
generalized Beth-Uhlenbeck form of the equation of state which is systematically extended from 
studying mesonic correlations [1] to the inclusion of diquark degrees of freedom [2, 3]. To that end 
we will employ a Nambu-Jona-Lasinio-type quark model with fourpoint interactions in mesonic 
(quark-antiquark) and diquark (quark-quark) channels. We shall discuss here the importance of the 
interplay of the resonant states with the residual non-resonant ones in the continuum of scattering 
states. Due to the Levinson theorem both contributions have the tendency to compensate each other 
in quark matter above the Mott transition [4], see also [2, 3, 5, 6]. 

The most intriguing questions will occur when on the basis of this in-medium bosonized ef¬ 
fective chiral quark model the next step of the hadronization program nucleon will be performed 
and diquarks will be “integrated out” in favor of baryons so that nuclear matter can be described 
in the model QCD phase diagram. A possible scheme for introducing baryons as quark-diquark 
bound states and integrating out the colored and therefore not asymptotically observable diquark 
fields has been suggesfed by Cahill and collaborafors [7, 8, 9, 10]. If was afterwards elaborafed by 
Reinhardf [11] and developed furfher by including fhe solifonic aspecfs of a field theoretic descrip¬ 
tion of the nucleon [12]. However, the step to describe nucleonic many-body systems on this quark 
model basis has not been performed by these authors. Exploratory studies within the framework 
of an effective local NJL-type model for the quark-diquark interaction vertex have revealed a first 
glimpse at the modification of the nucleon spectral function in the different regions of the model 
QCD phase diagram, including chirally restored and color superconducting phases [13]. Again, 
the crucial step towards an equation of state for nuclear matter with this microscopic model for 
nucleons in the QCD phase diagram has not yet been done. 

We will prepare the ground for a Beth-Uhlenbeck description of nuclear matter, to be discussed 
in future work. In particular, at zero temperature the structure of a Walecka model for nuclear matter 
shall emerge under specified conditions. Earlier work in fhis direcfion [14,15,16] has demonsfrafed 
fhis possibilify allhough no unified descriplion of fhe nuclear-lo-quark mailer Iransilion has been 
possible and fhe elucidalion of physical mechanisms for fhe very Iransilion belween fhe hadronic 
and fhe quark mailer phases of low energy QCD has been spared oul. Or sludy aims al indicating 
directions for filling fhis gap by extending fhe discussion of fhe Moll mechanism for fhe dissociation 
of hadronic bound slates of quarks wilhin fhe NJE model descriplion of low-energy QCD which has 
been given on the example ot two-particle correlations (mesons and diquarks) in [2] and sketching 
steps towards a Beth-Uhlenbeck description of nuclear matter where nucleons are treated as quark- 
diquark correlations in quark matter [17]. 

1.1 Quark matter thermodynamics in the NJL model 

We base the approach on the NJL model Lagrangian including vector and diquark interaction 
channels besides the standard scalar-pseudoscalar chirally symmetric interaction for the isospin 
symmetric case (/r„ = iXd = fi and = mf) 
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^ = (7[i^-mo + 7oAt]<? + -^int , (1-1) 

^int = Gs [{qqf+ {q^Y5T:qf] +Gy{qiyfiqf + G^ Yj {q^y5'^2^Aq^){q^^Y5'^2^Aq) ■ 

A=2,5,l 


Starting from the Lagrangian, we perform the usual bosonization by means of Hubbard- 
Stratonovich transformations, thus integrating out the quark degrees of freedom to obtain a path 
integral representation of the partition function (and thus the thermodynamical potential Q.) in 
terms of composite fields, mesons (M = a, 71, (O^) and (anti-)diquarks (D,D = A^, A^, A = 2,5,7), 


iT 




-fdV 


T^+;r^ 

4 Gs 


<^1 I |AaI^ 
4 Gv 4 Gd 


| + jlndet{/t5 * [frA.ffljiAAAll} 


( 1 . 2 ) 


where the inverse propagator is decomposed into mean field and fluctuation, S ^ = Sq ^ — Z, with 

^-1 ^ f{hn+l^*)yo-y-p-m AMFi75'^2A2 
2 ^ AJ;^pi75T2A2 {izn-l^*)yo-y-p-m 

a + iysf-K -AAiysTii^A 

-A^i75T2AA a + i75f^-7r 

Here, Zn = (2n + l)7rr are the fermionic Matsubara frequencies and we have introduced the dy¬ 
namical quark mass m = mof- a and the effective chemical potential jJ.* = jJ. — CUmf; the vector 
meson has been redced to the meanfield value of its timelike component co^ = (cUmfjO). In the 
present work, we consider NJL case, but we will restrict ourselves to the normal phase without 
color superconductivity {Amf = 0). Moreover, since we want to describe baryons as quark-diquark 
bound states, we have to go beyond Gaussian order in the diquark fields and include an infite set of 
diagrams which can be resummed to define the quark-diquark Bethe-Salpeter kernel, see [17] and 
references therein. 

As a result, the thermodynamic potential Q. =—{T/V)\nf^ takes the form 





^tot — ^MF A ^res ; 


(1.5) 


where the meanfield part Hmf = ^cond + consists of a contribution from condensates 


^cond — 


^MF _ ^F 

4Gs 4Gv 


(1.6) 


and from dynamical quarks 




Hq = -INcNf I f^\E„ + Tln 


The correlation contributions have the form 


ax = ±y^Trln(r^Sx‘) ,X = M, D, D, B, 


(1.8) 
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where for fermions ((B)aryons) holds the plus sign and cx = 1 while for bosons ((M)esons, (D)iquarks 
and their antiparticles (D)) holds the minus sign and cx = 2. The coeffhcient dx in (1.8) is the de¬ 
gree of degeneracy of the state X for which the inverse propagator can be given the form [2] 

*5x' = 7!--nx(iz«,q) (1.9) 

Dx 

with the effective coupling constant Gx being defined by the coupling in the corresponding in¬ 
teraction channel of the NJL Lagrangian. The coupling for the baryon Gb is discussed below. 
The polarization functions nx(iz«,q) are defined in fhe RPA approximation as one-loop infegrals 
which involve combinafions of quark-quark, quark-anfiquark and quark-diquark propagators wifh 
fhe corresponding vertex functions in fhe meson-, diquark-, and nucleon channels, respecfively. 
Any residual terms which cannot be captured in the above approximation scheme are summarized 
in the residual thermodynamic potential fires- Here one finds in particular confribufions beyond fhe 
Gaussian approximation in fhe meson fields as well as confribufions fo fhe baryon-baryon inferac- 
fion. The required polarizafion loop infegrals for mesons and diquarks are given in Ref. [2] for fhe 
NJL model. For fhe Polyakov loop extended NJL model one has fo replace in fhe final expressions 
fhe Fermi- and Bose disfribufion funclions for fhe color-carrying degrees of freedom (Q, D, D) by 
fhe generalized Fermi and Bose disfribufion functions, respecfively [3, 17, 18]. 


1.2 The baryon as quark-diquark state 

The quark-diqark polarization loop, which defines the baryon as an effective two-particle cor¬ 
relation in the medium, has in general a Dirac structure which, however, for a baryon at rest reduces 
to the scalar function (see also [17]) 


nB(iz„,0) = 4m 


d^p 


{Inf AEqEo 
/+(£,)+g-(£z)) 
iZn-E„+Eu 


\-r{Ef+g-{ED) 
T Eq “h Ed 

f-{Eq)+g+{ED) \ 

iZn+Eq-Eo J ’ 


1-/+(£,)+g+(£z)) 
iZn -Eq-E d 


( 1 . 10 ) 


where Eq = sj and Ed = y p^ + m^. For the PNJL model, the distribution functions for 
fermions {f'^{E) = l/[exp(£' =F jU)/r -|- 1]), bosons (g"'"(£') = l/[exp(£' =F/i)/r — 1]) and their 
antiparticles need to be replaced by their generalizations for the PNJL case and get then strongly 
suppressed in the confining phase. 



Figure 1: Diagrammatic representation of the instantaneous approximation to the quark-diquark interaction. 
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In order to find the baryon mass, we make the instantaneous approximation for the quark- 
diquark interaction vertex (see Fig. 1) which corresponds to identifying 


Gb 


^Sl-qq 

m 



Amo 

anD((3,o) I 

dti) \(i)—mD 


const, 


( 1 . 11 ) 


and solve the Bethe-Salpeter equation for the baryon as a bound quark-diquark state with a mass 
ms according to 


l-GBne(mB,0)=0. (1.12) 

Note that in contrast to Ref. [13] where Gb = const was assumed, we follow here [17] and keep 
the dependence on the quark mass which stems from the instantaneous limit of the quark exchange 
kernel, see Fig. 1 . Since with increasing temperature towards the chiral restoration the quark mass 
drops, this coupling gets enhanced and therefore the tendency for the nucleon being a “borromean 
state” which was already noted in [13] gets enhanced. 

We can solve the mean field gap equation for the quark mass as a function of temperature and 
chemical potential for the NJL model, which will serve as inputs for solving the Bethe-Salpeter 
equations for mesons, diquarks and baryons in the medium. The masses of these states are ob¬ 
tained as poles of their propagators and are shown in Fig. 2. We observe that the chiral symmetry 




Figure 2: Left panel: Mass spectrum of pions, sigma mesons, diquarks and baryons as functions of the tem¬ 
perature T. Also shown are the relevant thresholds: Iniq for mesons and diquarks and niq + nio for baryons. 
Right panel: Binding energy (mass defect) for baryons relative to the three-quark or the quark-diquark con¬ 
tinuum, resp., and for diquarks relative to the two-quark continuum as a function of the temperature T. The 
baryon in hot quark matter is a “borromean” state: when the diquark becomes unbound, the baryon binding 
energy is still nonvanishing. 

restoration which is reflected in the dropping quark mass function, induces a Mott effect for the 
pion and the scalar diquark. For both states the kernel of their Bethe-Salpeter equation contains a 
Pauli blocking term since they are composed of two fermions. This Pauli blocking partly compen¬ 
sates the effect of the quark selfenergies (dropping quark masses) and leads to a stabilization of the 
bound state masses against medium effects. This results in the crossing of the bound state masses 
with the continuum threshold, leading to the dissociation of these bound states. 
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Interestingly, the situation is quite different for the baryons. Here, the bound state is composed 
of a fermion and a boson, so that there is a partial compensation of the phase space occupation 
effects in the Bethe- Salpeter kernel: the Pauli blocking for quarks gets partly compensated by 
a Bose enhancement for the diquarks. As a consequence of this, the baryon mass follows the 
behaviour of the quark-diquark continuum and does not experience a Mott effect in the range of 
temperatures we studied here, see Fig. 2. In the right panel of Fig. 2, we show the mass defect 
(binding energy) for diquarks and baryons, where for the latter we show two alternative definitions, 
i.e., relative to the mass of three quarks and relative to the sum of the quark and the diquark mass. 
As a remarkable fact we obtain that in a temperature region where the diquark became already 
unbound due to the Mott effect, the baryon is still bound. One can therefore say that the baryon 
in the vicinity of the chiral restoration transition behaves like a “borromean state”: while the two- 
particle state (diquark) got dissociated, the three-particle state (baryon) is still bound. 

1.3 Thermodynamics of meson, diquark and baryon correlations in quark matter 

The derivation of the diquark thermodynamic potential follows Ref. [2] where it was given 
for the NJL model and the form of a generalized Beth-Uhlenbeck EoS was derived. In such a 
formulation, the dissociation of bound states in a hot, dense medium by the Mott effect is given by 
the behaviour of the in-medium phase shifts which encode the analytic properties of the propagator 
in the complex energy plane in a polar representation. The analoguous result is given in the Ref. [3] 
for the PNJL model. 

Now one is in the state to study the thermodynamics of the meson-, diquark-, and baryon 
correlations in a hot and dense medium. The equations of state of interest can be derived from 
the thermodynamical potential (1.5), e.g., for the pressure P = —fltot or the energy density e = 
fltot — TdQ.ioi/dT — whereby all corelation contributions have the form of a Beth- 

Uhlenbeck equation 


fix 



0 


(1.13) 


where the upper (lower) sign holds for fermions (bosons) and T± (ft)) = (l ± The 

corresponding degeneracy factors are = do = dj^ = 3 and ds = 4. Note that in evaluating (1.7), 
(1.13) the zero-point energy terms will be dropped ("no sea" approximation). 

For the phase shifts a general decomposition into a resonant (R) and a continuum (c) part can 
be made 


5x(ft),q) = 5x,R(ft),q) + 5x,c(ftJ,q) , (1-14) 

where both parts are uniquely defined by the propagator of the correlation, see [2]. The bound state 
mass is located at the jump of the phase shift from 0 to tt and this jump corresponds to a delta- 
function in the Beth-Uhlenbeck formulas for the correlations. In the case when the continuum of the 
scattering states is separated by a sufficient energy gap from the bound state, it can be neglected and 
we obtain as a limiting case the thermodynamics of a statistical ensemble of on-shell correlations. 
We shall make use of this fact in the next section when deducing the Walecka model. 
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2. The Walecka model as a low temperature limit of the NJL model 

In this section, we show that the Walecka model in mean field approximation is contained in 
our model as limiting case. It is well known, that quarks undergo a Cooper instability leading to di¬ 
quark condensates in this model due to the attractive quark-quark interaction. At zero temperature, 
this happens as /i* > m or, if diquarks are bound, for /r* > mull. Here, we restrict our consider¬ 
ations to the region fL* < min{m,m/)/2}. However, by construction, the baryon appears as bound 
state, meaning Mb/3 < min{m,m£)/2}. Consequently, there is a range in chemical potentials, 
where the appearance of baryons already affects the thermodynamics and leads to non-vanishing 
baryon densities, 

Mb/3 </I* < min{m,mz)/2} , (2.1) 

but quarks and diquarks are not excited yet. One could easily generalize this to finite temperatures, 
but then the effects of mesons should be taken into account. 


2.1 Deducing the equations 

Let us discuss this in more details. Our thermodynamic potential is given by 

^tot = ^cond + + ^M,R + ^M,c + ^D.R + ^D,c + ^D.R “h ^D,c “h ^B,R + ^B,c ; (2.2) 

with the explicit form of each term given above, and where we have decomposed the meson, 
diquark and baryon contributions in the respective resonance and continuum parts, as discussed 
above. At zero temperature, there are no contributions from the scattering states, due to the re¬ 
striction n* < m, which happen to be the thresholds defined in fhe imaginary parfs of fhe meson 
and diquark polarizafion funcfions. The same holds for fhe scattering par! of fhe baryons. This 
drasfically simplifies fhe calculation as many ferms drop ouf. However, fhe low femperafure region 
is already of inferesf for fhe properfies of neufron sfars. 

We choose fhe diquark coupling sfrengfh Gd small enough fo avoid fhe diquark pole, and 
Iherefore avoid colored objecfs in Ibis regime, fhe regime where hadrons are fhe only relevanf 
degrees of freedom. Consequenfly fhere is no bound sfafe confribufion from fhaf parf eifher. Thus, 
(2.2) reduces fo 


^tot — ^cond + + ^M.R + ^B,R • 


(2.3) 


For fhe quarks, only fhe vacuum part confribufes in confrasl fo baryons, due fo fhe condition 
Mb/3 < H*. Explicifly we have 

A 


f2tnr — 


CO, 


'0 


4Gs 4Gv 


- 2NrN 


ci^f 


(271) 


E„-2N^ 


3^p 


0 0 
The infegrals can be evaluafed analyfically, so fhaf we can wrife 


Ab , Pf 


M2+P2 


. (2.4) 


^tot - 


COq 1 


4Gs 4Gv 4712 


-^[3F(A,m)+F(AB,MB) + G(Pf,MB)] , 


(2.5) 


wifh funcfions F and G given in fhe appendix. Pf = y (Pb)^ ~^b Fermi momenfum of fhe 
baryon, /ig = 3/i*. 
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2.2 No sea approximation 

The transition to the original Walecka model needs a few more steps. First, we drop the 
vacuum contributions to the thermodynamic potential (2.5); this is what in the Walecka model is 
referred to as the no-sea approximation. The resulting thermodynamic potential reads 


f^tot - 



4Gv 


1 


G(Pf ,Mb) . 


( 2 . 6 ) 


This already looks like a Walecka type model, but the scalar and vector fields here are connected 
to the quarks. We now need to clarify the connection between the quark order parameters and the 
corresponding baryonic ones used in the Walecka model. This is easily achieved by employing 
quark counting rules, where the vacuum contributions are absorbed into the vacuum nucleon mass: 
Mb = 3(mo + a) — Kb = Mb, vac — 3amed- Motivated by our finding of Subsection 1.2 (see also 
Fig. 3), that the baryon mass defect (binding energy) in the model as a quark-diquark bound state is 
in lowest order medium independent, we adopt here in a first step a constant value of the parameter 
Kb = 261.7 MeV. The solution for the density dependent effective baryon mass is shown in Fig. 3. 
The expression for amed proportional to the scalar density is then obtained from minimizing the 



Figure 3: Baryon mass Mb and three-quark continuum (3m) as function of the baryon density ny. 


thermodynamic potential (2.6) with respect to the order parameter a and setting A = Ab, 


COo 

ny 

ns 

turned 


T-Gyny 
2 1 

1 . ^ Mb 

^F,(Pp,Mb) = ^ 

3Gs 


9 Pp PPp 

PpEp-Mlln 


Mb 


271^ 


p2(PF,MB) = 3(2Gs)ns 


(2.7) 

( 2 . 8 ) 

(2.9) 

( 2 . 10 ) 


We can rewrite the pressure P = —fltot and energy density £ in terms of densities, keeping in mind 
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that the baryon Fermi momentum is Pf = , 

P = Gynl-9Gsnl + ^(^PFEF (^-Pj-M^ (2.11) 

e = -P + pny (2.12) 

= Gyn\ + 9Gsnl + ^{pFEF {iPj+Ml) . (2.13) 

Note the factor 9, which appears in front of the scalar coupling strength. This reminds of the quark 
origin of this parameter and is in contrast to the vector coupling, where there is no prefactor, as the 
latter is entirely due to baryons. We use NJL parameters which would correspond to a constituent 




Figure 4: Comparison of the pressure and the binding energy per particle in a fitted Walecka model and an 
unmodified NJL parametrization (Gy = 4.5 Gs)- For the empirical value of the saturation density we take 
no = 0.16 fm^^. 


quark mass of m = 400 MeV [19], i.e., mo = 5.588 MeV, A = 587.9 MeV, GsA^ = 2.442 and 
Gy = 0.5Gs. Comparing this to the usual expressions for the Walecka model, we can define NJL 
model coupling strengths which would reproduce the Walecka model 


9Gs,Walecka 


1 gl ^ 

) Gy .Walecka 


2 ml ' 


(2.14) 


For values, reproducing the binding energy and saturation density [20]: ma = 550 MeV, mo, = 783 
MeV, g(j = 10.3 and = 12.7, we can compare with the NJL model parametrization and find 
instead 


Gs = 0.33 Gs.Walecka , Gy = 0.024 Gy.Walecka ■ (2.15) 

This leads already to saturating matter, even though at too high densities (« 3.4 nf) and with a too 
large binding energy (« —40 MeV instead of —16 MeV), see Fig. 4. 


3. Discussion and outlook 

The general aim of our work is the derivation of the thermodynamics for a meson-baryon 
system (quantum hadrodynamics = QHD) on the basis of quark degrees of freedom and to explore 
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the application limits of QHD, e.g., in the case of relativistic mean-field models for nuclear matter 
in heavy-ion collisions and neutron stars. The present contribution should be seen as a first step in 
this direction. 

We could show that, when the analysis is restricted to a mean-field approximation for mesons 
and to small densities below the quark threshold (/Tb < 3m^), the structure of the well-known a — CO 
model (Walecka model) for nuclear matter emerges in the limit T —> 0. In this way we obtain from 
the NJL model an equation of state for nuclear matter at finite densities and temperatures, which 
has a saturation property. Taking the model as it is and using the standard model parameters, 
however, we could not reproduce the phenomenological values for the saturation properties. In 
fact, taking into account the simplicity of the model, this was not to be expected. Nevertheless, it 
should be viewed as an achievement of this work that a general path from a chiral quark model to 
the thermodynamics of nuclear matter could be identified. 

In the present case, we find that the saturation density and the binding energy of the NJL nu¬ 
clear matter, which result from the interplay between the scalar and vector mean fields, are too 
large. In order to obtain the phenomenological values, one should take in to account additional ef¬ 
fects which lead to an enhancement of both, the repulsion (vector field) and of the attraction (scalar 
field). Here we think that the consideration of higher-order terms in the meson fields could be a 
promising option, because a consistent bosonization of the NJL model requires the consideration 
and renormalization of terms up to the fourth order, which would otherwise lead to divergencies. 
In this way, one obtains, for example, from the NJL model the successful hadronic linear a model, 
see [21] and references therein. Alternatively, terms of higher order can already be introduced on 
the level of the original quark Lagrangian, leading, however, to new coupling constants (see [22] 
and refs, therein). This freedom has been used, e.g., by Huguet et al. [23] to successfully reproduce 
the saturation properties of nuclear matter at T = 0 in this way. Other interesting approaches to 
describe nuclear matter based on quark degrees of freedom are discussed in [2]. Besides this, we 
have omitted terms in Ores which correpond to quark exchange terms between baryons and would 
result in a repulsive Pauli blocking contribution at high densities. Considering such higher order 
terms one should eventually be able to reduce the binding energy and saturation density towards 
their physical values! 

Finally, we want to point out that a suitable confining mechanism has to be present in a quark- 
matter model, in order to exclude the appearance of unphysical degrees of freedom already in the 
experimentally well studied nuclear-matter regime. In this context the Polyakov loop has proven to 
be an interesting tool, although it does not involve all aspects of confinement. This is in particular 
the case at low temperatures and for T = 0. 
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A. Basic integrals 

We list some integrals encountered in the zero temperature evaluation of quark loop integrals 
F{x,y) = S J dp \J +3;2 = +y^{2x^ +y'^) -y^ln ^~^ ^ (A.l) 

Partial derivatives of these functions are encountered in the evaluation of generalized gap equations; 
here we require 


dF{x,y) 

F2[x,y) = —— =4y 


/ 2 I 2 2, x+y/x^+y^ 

xv-^ +y — y In-- 


(A.3) 
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